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ON  THE  ESTIMATION  OF  AREAL  WINDSPEED  DISTRIBUTION 
IN  TROPICAL  CYCLONES  WITH  THE  USE  OF  SATELLITE  DATA 


Andrew  Timchalk 
Meteorological  Satellite  Laboratory 
National  Environmental  Satellite  Service 
NOAA,  Washington,  D.C. 


ABSTRACT.  This  paper  Is  an  attempt  to  determine  whether 
the  areal  windspeed  distributions  In  tropical  cyclones 
can  be  estimated  by  use  of  two  satellite  data  parameters: 
(1)  the  low-level  convectlve  cloud  band  crossing  angle 
relative  to  tangents  to  concentric  circles  centered  on 
the  eye,  (2)  the  Infrared  temperatures  from  the 
10.5-to-12.5  \m   channel.  These  two  parameters  plus  the 
"ground  truth"  winds  were  all  derived  for  128  polar 
coordinate  grid  points  (16  azimuths  and  8  radii) . 

Coincident  values  of  the  variables  were  evaluated 
statistically  first  by  linear  regression  and  afterward 
by  multiple  screening  regression.  For  the  entire  storm 
area  (all  radii) ,  the  correlation  between  observed  wind- 
speed  and  the  two  satellite-derived  parameters  [ (1)  and 
(2),  above]  was  0.4  and  -0.6,  respectively.  However, 
when  evaluated  at  separate  radii  the  two  satellite 
parameters  were  poorly  correlated  to  either  the  observed 
or  relative  windspeed. 

Radial  distance  from  the  storm  center  was  found  to 
explain  63%  of  the  variance  In  windspeed.  When  this 
relation  was  removed  by  screening  regression,  the  two 
satellite  parameters  could  explain  only  an  additional 
2%  to  3%  of  the  variance. 

Correlation  between  the  cloud  band  crossing  angle  and 
the  wind  field  crossing  angle  was  less  than  0.5  at  all 
radii.  Correlation  between  the  windspeed  and  the  wind 
field  crossing  angle  was  less  than  0.2.  In  both  cases, 
the  coefficients  were  lower  than  expected. 

I.   INTRODUCTION 

A.  Problem  and  Scope 

The  purpose  of  this  work  is  to  determine  whether  the  areal  windspeed 
distributions  in  tropical  cyclones  can  be  estimated  using  satellite  data 
parameters.  The  scope  of  this  study  is  confined  to  two  such  parameters. 
Low-level  convectlve  cloud  band  crossing  angles  relative  to  tangents  to 


concentric  circles  centered  on  the  eye  were  obtained  from  the  0.5-to-0.7  pm 
visual  channel  and  infrared  brightness  temperatures  from  the  10.5-to-12.5 
ym  infrared  channel. 

Because  of  expected  vertical  wind  shear  at  radii  beyond  100  n.mi.  from  a 
storm  center,  aerial  reconnaissance  winds  are  probably  unrepresentative  of 
winds  near  the  surface.   Therefore,  this  study  was  designed  to  estimate  the 
distribution  of  the  maximum  sustained  windspeed  at  ship  deck  level. 
Satellite  parameters  and  wlndspeeds  were  derived  over  a  128-point  polar 
coordinate  grid  consisting  of  points  located  at  16  azimuth  angles  and  8 
radii  from  storm  center. 

B.  Past  Studies  of  Wind  Distribution 

Prior  to  1943,  when  the  military  services  of  the  United  States  began 
aerial  reconnaissance  into  tropical  cyclones,  practically  all  the  data  used 
for  the  study  of  wind  distributions  in  storms  were  measured  at  island  and 
coastal  stations.   Based  on  empirical  observations  indicating  that  angular 
momentum  was  only  partially  conserved,  Byers  (1944)  predicted  that  the  wind- 
speed  should  be  inversely  proportional  to  the  square  root  of  the  radial 
distance  from  the  storm  center.   Similar  wind  structure  in  tropical  cyclones 
has  since  been  noted  by  numerous  investigators. 

The  "classical"  model  of  the  wind  field  surrounding  tropical  cyclones  was 
developed  by  Hughes  (1952).   By  using  low-level  (about  1,000  feet)  recon- 
naissance data  for  the  period  1945-47,  Hughes  derived  a  typical  windspeed 
distribution  for  a  large,  stationary,  mature  cyclone.  Hughes  found  that 
when  he  substituted  the  tangential  wlndspeeds  at  radii  1°  and  3°  (60  and 
180  n.mi.)  into  the  equation: 

Vx  =  c  CD 

(where  R  is  the  radius,  v«  the  tangential  windspeed,  C  and  x  are  constants), 
x  =  0.62  provided  the  best  fit.   This  compares  with  the  0.5  predicted  by 
Byers  (1944)  and  Riehl  (1963) . 

Jordan  (1952)  constructed  a  picture  of  the  upper-level  wind  circulation 
around  tropical  cyclones.  Miller  (1958)  investigated  the  three-dimensional 
wind  structure  around  hurricanes.  Hawkins  (1962)  examined  the  vertical 
profiles  of  horizontal  wind  in  hurricanes  and  found  that  the  preservation 
of  the  windspeed  with  height  is  well-marked  below  20,000  feet  within  a  radius 
2.5  times  the  radius  of  maximum  winds.   These  and  similar  studies  have  not 
only  provided  data,  but  have  also  served  as  a  basis  for  many  subsequent 
theoretical  and  experimental  studies. 

A  search  of  the  literature  shows  that  most  studies  of  wind  distributions 
are  either  of  Individual  storm  cases  or  are  of  the  high  energy  core  of  the 
storms;  that  is,  within  radii  of  less  than  100  n.mi.  from  the  storm  center. 
This  limitation  is  mainly  due  to  the  lack  of  low-level  aircraft  recon- 
naissance beyond  the  100  n.mi.  (Shea  and  Gray  1972). 


C,  Need  for  the  Study 

Knowledge  of  the  low-level  wlndspeed  In  tropical  cyclones  Is  Important 
for  the  safety  of  all  ocean  going  vessels.   The  areal  distribution  of  wind- 
speed,  e.g.,  the  extent  of  the  30-  and  50-kt  isotachs  is  important  in 
preparedness  or  evacuation  decisions  for  coastal  installations  (Brand  and 
Blelloch  1975).  Distributions  of  wlndspeed  and  inflow  angle  affect  the 
height  of  storm  surges  that  accompany  intense  tropical  cyclones  (Jelesnlanski 
1966) ;  and  areal  speed  distributions  are  of  primary  importance  in  the  study 
of  the  dynamics  and  energetics  of  these  storms. 

Weather  satellite  capability  to  view  large,  remote  areas  of  the  tropics 
simultaneously,  and  the  increasing  cost  of  aircraft  reconnaissance,  are 
compelling  reasons  to  search  for  methods  of  specifying  the  wlndspeed  in 
tropical  cyclones. 

Satellite  data  are  being  used  effectively  to  locate,  and  to  estimate  the 
maximum  intensity  of,  tropical  cyclones  (Sheets  and  Grieman  1975);  however, 
little  study  has  been  devoted  to  the  use  of  these  data  in  specifying  the 
areal  distribution  of  the  wlndspeed  in  storms. 

II.  METHOD 

A.  Selection  of  Satellite  Data  Parameters 

1.  Selection  of  Cloud  Band  Crossing  Angle. 

A  hurricane  is  a  thermally  driven  circulation  whose  primary  energy  source 
is  the  release  of  latent  heat  of  condensation.  A  model  hurricane  shows  a 
wind  field  with  cyclonic  inflow  in  the  lower  troposphere  and  anticyclonic 
outflow  aloft.  Riehl  (1963)  obtained  some  reasonable  results  from  a  simple 
two-layer  model  he  developed  in  which  the  storm  is  assumed  to  be  symmetri- 
cal and  steady-state.   Conservation  of  absolute  angular  momentum  in  the 
outflow  layer  and  conservation  of  potential  vortlcity  in  the  inflow  layer 
were  the  key  assumptions. 

Through  the  use  of  the  absolute  angular  momentum  equation,  the  equation 
for  the  tangential  component  of  motion,  and  assuming  conservation  of 
potential  vortlcity  in  the  inflow  layer,  Riehl  derived  an  expression  for 
the  tangential  wlndspeed  distribution  in  tropical  cyclones  eq.  (11) .  His 
initial  formulation  contained  an  inflow  angle  term  (cos  a) ,  where  a  is  the 
angle  the  wind  makes  with  a  tangent  to  a  circle  centered  on  the  storm. 
This  term  was  later  omitted  after  applying  the  findings  of  Ausman  (1959) 
because  the  value  of  cos  a  was  near  unity  when  the  crossing  angles  were 
small  (i.e.,  nearly  tangential  flow).  However,  this  term  will  be  retained 
following  the  derivation  of  Riehl  (1963). 

The  equation  of  motion  in  polar  coordinates  for  the  tangential  compo- 
nent, neglecting  lateral  friction,  is: 


dv0    v0vr  1  3P     1  3Te, 

— 2-  +  -2-1  +  fv  .  -  +  22.        (2) 

dt     R      r     p  36     p  3z         u; 


where  v©  is  the  tangential  velocity,  vr  the  radial  velocity,  and  Tqz  the 
tangential  shearing  stress  in  the  0-z  plane.  The  absolute  angular  momentum 
per  unit  mass  about  a  vertical  axis  at  the  storm  center  is  given  by: 

£R2 
fl  -  veR  +  -—■  .  (3) 

Differentiating  eq.  (3)  with  respect  to  time  and  noting  that  dR/dt  =  vr, 
yields: 


dQ  dvo 

=  R — D 

dt     dt 


=  R  -r2-  +  vQvr  +  Rfvr  .  (4) 


Multiplying  eq.  (2)  by  R,  eq.  (4)  can  be  written  as: 

JL.  _  R  3p  ,  R  3l6z  . 

dt      p  36   p  3z 

Assuming  a  symmetric  vortex,  3P/36  «  0,  so 


(5) 


dJ3     R  3T 

■  mmmm         SZ       wmmmmmm 

dt     P  3z 


0z 


(6) 


Eq.  (6)  shows  that  the  time  rate  of  change  of  absolute  angular  momentum 
depends  on  the  vertical  gradient  of  the  shearing  stress. 

For  the  inflow  layer,  eq.  (6)  can  be  integrated  over  an  inflow  depth  of 
6z  where  at  the  top,  Tqz  =  0. 

/  — P<5z  "  -T0s  '  (7) 

at 

where  Tq8  Is  the  surface  tangential  shearing  stress.  Assuming  conservation 
of  potential  vorticity,  the  curl  of  the  component  of  frlctional  force  along 
6  must  be  zero.   That  is,  V  x  F«  =  0  where 

Fq  ■  fiz.  Thus  for  FgR  =  constant  (and  by  integrating  from  the  surface 

p  dz 


to  the  height  where  the  stress  vanishes)  the  result  is: 

RT0S  =  constant  .  (8) 

At  this  point  an  empirically  derived  expression  for  the  surface  stress  is 
introduced : 

Tes  =  cdpsve/cos  a  *  (9) 

Cj  is  the  drag  coefficient,  ps  the  surface  density,  vqs  the  surface  tan- 
gential windspeed  (at  ship  deck  level)  and  a  is  the  inflow  or  crossing 
angle. 

Combining  eq.  (8)  and  (9)  and  solving  for  V0g: 


)s 


1uos  cx/r;  ,  qq) 


v0^  =  C,  (cos  a/R)'5 
i 


where  C-^  =  (constant /C^ps)  .  Eq.  (10)  shows  that  the  tangential  windspeed 
is  directly  proportional  to  the  square  root  of  cos  a.  Hence,  windspeed  is 
inversely  related  to  the  crossing  angle.  Note  that  when  the  flow  is 
tangential,   i.e.^a  =  0,  cos  a  ■  1.  Thus  eq.  (10)  reduces  to  the  well- 
known  tangential  windspeed  distribution  equation  where  the  subscript  (s) 
has  been  omitted  for  convenience: 


h 


vj?   =  Cn  .  (ID 

For  convenience  also,  eq.  (11)  will  be  written  in  the  form: 

v6  =  CxlT  *5#  (12) 


In  hurricanes  the  tangential  component  of  the  windspeed  (v  )  is  usually 
much  larger  than  the  radial  component  (vr) ;  therefore  the  total  windspeed 
should  be  closely  approximated  by  the  tangential  component.  Further,  if 
the  cloud  bands  are  alined  with  the  wind  direction,  then  as  indicated  by 
eq.  (10),  the  observed  wind  speed  should  be  directly  proportional  to  the 
square  root  of  the  cosine  of  the  low  cloud  band  crossing  angle.  The 
relation  implied  by  eq.  (10)  along  with  satellite  studies  by  Timchalk  et  al. 
(1965),  Fritz  et  al.  (1966),  and  Dvorak  (1973)  were  the  bases  for  selec- 
ting the  crossing  angle  of  cloud  bands  (acid)  as  one  of  the  satellite  data 
parameters.  A  necessary  assumption  is  that  acid  is  nearly  equal  to  the 
crossing  angle  of  the  observed  wind  (otODS) ;  this  assumed  relation  will  also 
be  examined. 


2.  Selection  of  Infrared  Radiation  Temperature 

Riehl  (1954) ,  citing  high-level  observations  and  the  theoretical  work  of 
Haurwitz  (1935) ,  noted  that  mature  storms  extend  through  the  troposphere 
to  near  the  100-mb  level.  At  these  levels,  the  cloud-top  temperatures 
would  be  as  low  as  -78°C  (Schacht  1946).   Accordingly,  very  cold  cloud-top 
temperatures  should  be  found  over  intense  storms. 

The  work  of  Fritz  et  al.  (1966)  with  the  TIROS  visible  channel  data 
showed  a  positive  correlation  between  the  size  of  the  hurricane  "canopy" 
and  the  maximum  windspeed  in  tropical  cyclones.  Further,  in  an  attempt  to 
establish  a  relation  between  the  equivalent  blackbody  radiation  tempera- 
ture (Tug)  and  the  intensity  of  tropical  cyclones,  Hubert  et  al.  (1969) 
found  the  canopy  sizes  could  be  closely  approximated  by  the  extent  of  the 
260  K  infrared  isotherms.  They  noted  that  Tgg  at  the  wall  cloud  radius, 
where  the  maximum  windspeeds  are  located,  was  In  most  cases  colder  than 
230  K,  while  in  the  outer  areas,  where  the  windspeeds  are  low,  Tjjb  was  near 
300  K.  These  findings  suggested  that  windspeed  should  be  Inversely  related 
to  Tbb»  particularly  over  Intense  tropical  cyclones.  Because  satellite 
infrared  temperature  had  been  observed  to  be  related  to  the  windspeed,  TgB 
was  selected  as  a  study  parameter. 

B.  Analysis 

1.  Determination  of  Wind  Field 

Since  the  objective  of  this  study  was  to  estimate  the  areal  windspeed 
distribution,  it  was  necessary  to  obtain  the  "ground  truth"  wind  field. 
Following  the  recommendation  of  Frank  (1974) ,  aircraft  reconnaissance  data 
were  collected  via  the  Navy's  Hurricane  Circuit  located  at  the  Fleet 
Weather  Central  Suitland,  Md.  over  the  Atlantic  Ocean  during  the  1974 
tropical  storm  season.  Aircraft  reports  at  or  below  the  700-mb  level, 
rawinsonde  wind  data  for  the  609-m  level,  and  ship  and  buoy  data,  were 
used.  The  609-m  level  was  chosen  after  noting  the  work  of  Hawkins  (1962) 
and  examining  vertical  profiles  of  windspeed. 

All  the  data  were  composited  about  the  1600  GMT  position  of  the  storm. 
To  do  this,  the  bearing  and  distance  from  the  storm  center  of  each  obser- 
vation was  calculated  from  the  interpolated  storm  position.   The  obser- 
vations for  each  day  were  then  plotted  on  a  Mercator  projection  showing 
the  wind  direction,  speed,  type  and  height  of  wind  observation,  and  the 
time  difference  from  1600  GMT.   Only  wind  observations  that  were  made 
within  +  8  hours  of  1600  GMT  were  used. 

Applying  the  results  of  Riehl  (1954),  Miller  (1958) ,  Hawkins 
(1962),  and  others,  lsogon/ streamline  and  isotach  analyses  were  made, 
considering  such  factors  as  type,  height,  and  reliability  of  wind  observa- 
tion coupled  with  time  and  spatial  consistency.   These  field  analyses  were 
used  to  extract  the  wind  direction  and  speed  for  the  polar  coordinate  grid 
at  16  azimuth  angles  and  radial  increments  of  1*  latitude  (60  n.mi.)  out- 
ward to  a  distance  of  8°  latitude  (480  n.mi.)  from  "the  storm  center. 


Figure  1  shows  the  lsogon-streamllne  analysis  and  figure  2  shows  the 
Isotach  analysis  for  hurricane  Carmen  on  Sept,  1,  1974.  From  these  anal- 
yses, wind  direction  and  speed  were  extracted  at  each  of  the  128  grid 
points  shown  In  figure  3.  Wind  fields  were  prepared  for  10  days  selected 
from  3  hurricanes.   Copies  of  the  wind  field  may  he  obtained  by  contacting 
the  author. 

Since  none  of  the  10  storm  cases  was  stationary  over  the  16-hour 
composite  period,  smoothed  track  positions  were  used  to  calculate  a  mean 
storm  motion  vector  which  Is  printed  at  the  center  of  the  polar  grid  for 
each  case.  The  first  three  digits  indicate  the  azimuth  toward  which  the 
storm  was  moving,  while  the  remaining  two  digits  indicate  the  speed  of 
translation  in  knots.  A  list  of  the  storm  cases,  dates,  positions,  and 
storm  motion  vectors  used  in  this  study  is  given  in  table  1. 


Table  1.—  List  of  storm  cases 


Storm 

motion  vector 

Hurricane 

Date 

Location  at  1600  GMT 

direction/speed (kt) 

Carmen 

9/01/74 

17. 6N 

83. 9W 

274/13 

it 

9/02/74 

18. 8N 

88. 7W 

280/10 

ii 

9/05/74 

21. 5N 

90. 7W 

010/04 

ii 

9/06/74 

23. 5N 

90. 4W 

008/09 

ii 

9/07/74 

27. 5N 

90. 4W 

360/10 

ii 

9/08/74 

30, 3N 

92. 5W 

320/10 

Flfl 

9/17/74 

16. 9N 

80. 8W 

260/08 

it 

9/18/74 

16. 3N 

84. 3W 

270/15 

ii 

9/19/74 

16.  IN 

88. 0W 

290/08 

Camille 

8/17/69 

27. 5N 

88. 5N 

344/12 

In  moving  storms,  the  wind  at  a  given  point  contains  both  rotational 
and  translatlonal  components.   If  the  translation  component  (the  storm 
motion  vector)  Is  subtracted  vectorially  from  the  wind  at  each  grid 
point,  the  new  field  is  called  the  "relative"  wind  field.  Both  the 
observed  and  relative  wind  fields  were  studied  for  possible  relations 
with  the  two  satellite  parameters.  An  example  of  the  relative  wind  field 
for  hurricane  Carmen  on  Sept.  1,  1974  is  shown  in  figure  4.  Comparison 
of  figures  3  and  4  shows  the  difference  between  the  relative  and  the 
observed  wind  fields  for  a  storm  moving  toward  the  west  at  13  kt. 


SSE- 


Figure  1. — Isogon-streamllne  analysis  for  hurricane  Carmen,  1600 
Sept.  1,  1974.   Eye  position:   17. 6N,  83. 9W. 


GMT, 


I  "1^ 


Figure  2. — Isotach  analysis  (kt)  for  hurricane  Carmen,  1600  GMT,  Sept.  1, 
1974.  Eye  position:  17. 6N,  83. 9W. 
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Figure  3. — Observed  wind  field  (dddff)  of  hurricane  Carmen,  1600  GMT, 
Sept.  1,  1974;  windspeed  in  knots.  The  central  value  is  the  storm 
motion  vector. 
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Figure  4.— Relative  wind  field  (dddff)  of  hurricane  Carmen,  1600  GMT, 
Sept.  1,  1974;  windspeed  in  knots. 
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2.  Measurement  of  Cloud  Band  Crossing  Angle 

Except  for  the  hurricane  Camllle  case,  pictures  from  the  visible  channel 
CO. 5  to  0.7  um)  of  the  Synchronous  Meteorological  Satellite  (SMS-1)  were 
used  to  determine  the  orientation  of  the  low-level  convectlve  cloud  bands. 
The  resolution  of  the  SMS-1  pictures  used  here  was  1  n.ml.   For  Camllle 
ESSA  9  data  were  used,  with  a  picture  resolution  of  2  n.ml.  which  was 
quite  adequate  for  determining  the  orientation  of  the  cloud  bands. 
Figure  5  is  the  SMS-1  picture  of  hurricane  Carmen  on  Sept.  1,  1974,  while 
figure  6  shows  the  orientation  of  the  cloud  bands. 

Once  the  cloud  band  orientations  were  transferred  from  the  gridded  pic- 
tures to  a  Mercator  map,  a  polar  grid  overlay  was  used  to  extract  otcu 
at  the  grid  points. 

Unlike  a  wind  vector,  the  orientation  of  a  cloud  band  can  be  assigned 
two  directions,  as  shown  in  figure  7.   In  order  to  preserve  inward  (-) 
and  outward  (+)  measurements  of  the  crossing  angles,  the  direction  was 
always  assigned  in  a  cyclonic  sense  about  the  storm  center.  The  angles 
were  measured  with  respect  to  tangents  to  concentric  circles  centered  on 
the  storm  eye.  Note  that  the  limits  of  aci<j  are  +  90*.  For  points  where 
no  low  cloud  band  could  be  seen,  or  the  interpolated  distance  between 
cloud  bands  was  greater  than  3°  latitude  (180  n.ml.),  or  the  extrapolated 
distance  from  the  nearest  cloud  band  was  greater  than  1°  latitude 
(60  n.mi.),  999  was  encoded  to  indicate  that  the  crossing  angle  was  missing. 
Values  for  &cid   were  obtained  at  77%  of  the  polar  grid  points.  Missing 
values  were  usually  located  at  the  outer  radii,  over  land,  and  in  clear 
to  partly-cloudy  areas. 

3.  Extraction  of  Infrared  Radiation  Temperature 

Scanning  Radiometer  Infrared  (SRIR)  measurements  of  equivalent  black- 
body  temperature  (Tgg)  obtained  from  polar-orbiting  satellites  are 
routinely  mapped  in  digital  form  at  full  resolution  (approx.  5  n.mi.  at 
subpoint)  and  stored  on  magnetic  tape.  No  SRIR  data  tape  existed  for 
Camllle,  while  the  data  tape  for  Carmen  on  Sept.  5  contained  much  noise 
and  missing  data.   Thus  of  the  10  storm  cases,  SRIR  coverage  was  avail- 
able for  only  8.  Because  the  infrared  temperatures  (radiances)  were 
obtained  by  the  NOAA-II  polar-orbiting  satellite,  the  measurements  were 
made  at  slightly  different  times.  However,  all  the  measurements  were 
made  within  +  1  hour  of  the  1600  GMT  compositing  time  with  the  exception 
of  the  Sept.  1,  hurricane  Carmen  case  in  which  the  measurements  were 
made  at  1442  GMT.  This  case  is  shown  in  figure  8. 

A  program  was  devised  to  calculate  the  latitude  and  longitude  of  each 
grid  point;  Tgg  was  then  extracted  from  the  tape  for  these  locations. 
To  insure  that  the  Tgg  selected  would  represent  a  suitable  area,  Tjjb 
was  averaged  over  3x3  scan  spots,  or  a  15  x  15-n.ml.  area  for  each 

grid  point. 
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Figure  5. — SMS-1  picture  of  hurricane  Carmen,  1600  GMT,  Sept.  1,  1974. 
Eye  position:  17. 6N,  83. 9W. 
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Figure  6, — SMS-1  picture  of  hurricane  Carmen,  1600  GMT,  Sept.  1,  1974 

with  orientations  of  low-level  convective  cloud  bands  used  to  determine 
the  cloud  band  crossing  angle. 
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Axis  of  cloud  band  with 
crossing  angle  (ctcld)<0< 


Axis  of  cloud  band  with 
crossing  angle  (0c|cf)>o< 

Note:    +90°^ac|c|^_90o 


Figure  7. — Schematic  showing  how  the  cloud  band  crossing  angles  (acld) 
were  assigned. 
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Figure  8. — NOAA  II  Infrared  temperatures  (TBB)°K  for  hurricane 
Carmen,  1442  GMT,  Sept.  1,  1974. 
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Comparison  with  full-scale  resolution  data  indicated  that  the  3x3 
scan-spot  average  preserved  the  variations  in  temperature  over  the 
effective  radiating  surfaces,  and  was  judged  to  be  compatible  with  the 
scale  and  reliability  of  the  windspeed  data. 

C.   Correlation  Techniques 

Linear  correlation  and  regression  coefficients  were  calculated  between 
the  two  satellite  data  parameters  and  (first)  the  observed  windspeed  (vQ^s), 
and  (second)  the  relative  windspeed  (vrei) .   Correlation  coefficients  were 
calculated  separately  for  each  storm  case  at  each  radius.   Radially  cum- 
ulative correlations  for  each  storm  case  were  also  calculated  for  radii 
1°  to  8°.   Finally,  radial  and  radially  cumulative  correlations  were 
calculated  for  all  the  storm  cases  combined. 

Linear  correlations  were  also  made  between  the  observed  and  relative 
windspeeds  (v0bSj  vrel)  anc^  their  respective  crossing  angles  (a0bs,  a^ei)  • 

All  the  linear  correlations  were  evaluated  to  determine  whether  the 
parameters  were  related  to  the  windspeed  when  the  data  were  stratified 
either  by  storm  case  or  by  radius.   Both  satellite  parameters  were  found 
to  be  related  to  the  windspeed  and  to  the  radial  distance  from  the  storm 
center  (R) .   To  isolate  the  effectiveness  of  the  two  satellite  data 
parameters  in  specifying  the  windspeed  distribution,  a  Stagewise  Screening 
Regression  Method  and  later  a  Stepwise  Multiple  Regression  Method  were 
applied  to  the  data.   Results  of  the  correlation  and  regression  analysis 
are  discussed  in  section  III. 

D.   Data  Accuracy  and  Representativeness  of  Cases 

Mean  errors  were  estimated  to  be  +  5°  for  ac±d  and  +  3°C  for  the  full 
resolution  Tgg.  Natural  variations  of  the  wind  in  time  and  space  along 
with  inaccuracies  in  wind  measurement  make  it  difficult  to  assign  limits 
of  accuracy  to  the  observed  wind  fields.  Noting  the  natural  variations 
found  by  Gentry  (1964),  mean  errors  of  +  10°  in  direction  and  +  20%  in 
speed  may  be  too  small.  On  the  other  hand,  it  is  likely  that  the  large 
variability  found  by  Gentry  (ibid.)  has  been  damped  by  smoothing  in  the 
isogon-isotach  analyses. 

In  the  10  storm  cases  studied  here,  the  maximum  observed  surface  wind- 
speed  varied  from  70  kt  in  hurricane  Carmen  (9/05/74)  to  140  kt  in 
hurricane  Camille  (8/17/69).   Because  all  the  storm  cases  were  either  in 
the  Gulf  of  Mexico  or  the  Caribbean  Sea,  some  part  of  the  storm  was 
nearly  always  over  land.   However,  because  of  the  way  in  which  the  ob- 
served winds  were  composited,  it  is  unlikely  that  land  effects  distorted 
the  wind  fields  appreciably.   Further,  since  ac±^   was  usually  missing 
over  land  masses,  this  parameter  was  probably  not  significantly  affected 
by  the  proximity  of  land. 

In  an  attempt  to  isolate  any  nonrepresentative  cases,  linear  correla- 
tions were  first  calculated  by  individual  storm  case.   Although  there 
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were  marked  differences  between  cases,  no  pattern  of  differences  was 
detected,  possibly  because  of  the  small  number  of  cases.  Accordingly,  the 
individual  cases  will  be  discussed,  but  the  statistical  results  will  not 
be  presented  individually,  except  for  the  hurricane  Carmen  case  of  Sept. 
1,  1974,  which  will  be  presented  for  illustration. 

III.   RESULTS 

A.  Linear  Regression  of  Windspeed  on  Cloud  Band  Crossing  Angle 

1.   Observed  Windspeed  vs.  Cloud  Band  Crossing  Angle 

Table  2  shows  the  results  of  the  linear  regression  of  v0DS  on  algebraic 
values  of  otci(j  for  each  radius  (R)  for  the  Sept.  1,  1974  hurricane  Carmen 
case.   In  testing  the  relation  involving  the  crossing  angles,  correlations 
were  made  by  treating  the  angles  first  as  algebraic  and  then  as  absolute 
values.  The  algebraic  means  of  OLcid   became  more  inward  while  the  mean 
windspeed  decreased  with  increasing  radius  in  the  storm.   (Winds  were  used 
only  if  coincident  with  ol.-^  values.) 

Table  2. — Linear  regression  of  observed  windspeed  (v0^s)  on  algebraic 
values  of  the  cloud  band  crossing  angle  Cotc^<j)  at  separate  radii  for 
hurricane  Carmen,  1600  GMT,  Sept.  1,  1974. 


n 

R  (d< 

^g) 

v0bs  (kt> 

acld  (deg) 

r 

a  (kt) 

b  (kt/deg) 

S  (kt) 

16 

1 

43.9 

-17.81 

-0.40 

36.08 

-0.44 

10.87 

16 

2 

25.8 

-22.19 

-0.38 

18.00 

-0.35 

8.26 

16 

3 

20.4 

-26.56 

-0.16 

17.76 

-0.10 

6.56 

9 

4 

17.8 

-36.39 

0.45 

32.81 

0.41 

6.27 

8 

5 

16.4 

-42.50 

0.65 

30.31 

0.33 

4.08 

8 

6 

15.7 

-38.13 

0.08 

16.78 

0.03 

4.34 

7 

7 

14.4 

-35.71 

-0.64 

9.59 

-0.14 

2.68 

U 

8 

13.0 

-50.63 

-0.77 

0.14 

-0.24 

2.36 

The  correlation  coefficients  for  individual  radii  should  be  considered 
only  as  crude  estimates  because  the  sample  is  so  small.   For  a  given 
radius,  the  regression  equation  is:  v0fjS  =  a  +  bta^j];  the  standard 
error  about  this  line  is  given  under  S. 

Since  acid  for  the  Sept.  1  hurricane  Carmen  case  showed  only  inward 
cloud  band  crossing  angles,  the  linear  regression  coefficients  derived 
using  the  absolute  values  of  &q\&   are  not  shown  because  the  results  would 
be  identical  to  those  given  in  table  2.   In  cases  where  some  outward  -tadd 
values  were  measured,  the  regression  results  obtained  by  using  the  abso- 
lute values  would  be  slightly  different.   The  correlations  of  v0^s  with 
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CX  ,  ,  at  individual  radii  for  the  other  nine  cases  were  small  and  even 

C-Ld 

changed  sign  from  one  case  to  another. 

Table  3  shows  the  results  of  the  radially  cumulative  regression,  from 
R  ■  1°  to  8°,  of  v0|jS  on  the  algebraic  value  of  cxcid  for  the  Sept.  1, 
hurricane  Carmen  case.   The  bottom  row  shows  the  results  for  the  entire 
storm  area. 


Table  3. — Radially  cumulative  linear  regression  of  observed  windspeed 
(v qIjq)  on  algebraic  values  of  the  cloud  band  crossing  angle  (aci(j) 
for  hurricane  Carmen,  1600  GMT,  Sept.  1,  1974. 


n 

£R     (deg) 

vQbs   (kt) 

<*cld   (de8> 

r 
-0.40 

a   (kt) 

b   (kt/deg) 

S    (kt) 

16 

1-1 

43.9 

-17.81 

36.08 

-0.44 

10.87 

32 

1-2 

34.8 

-20.00 

-0.16 

30.73 

-0.20 

13.69 

48 

1-3 

30.0 

-22.19 

0.01 

30.37 

0.02 

13.76 

57 

1-4 

28.1 

-24.43 

0.18 

33.07 

0.20 

13.46 

65 

1-5 

26.6 

-26.65 

0.30 

34.78 

0.31 

12.90 

73 

1-6 

25.5 

-27.91 

0.33 

34.51 

0.32 

12.52 

80 

1-7 

24.5 

-28.59 

0.32 

32.92 

0.30 

12.43 

84 

1-8 

23.9 

-29.64 

0.34 

32.93 

0.30 

12.26 

Note  that  the  correlation  coefficients  change  sign  and  remain  positive  as 
the  radius  increases.   This  change  in  sign  appeared  in  all  but  one  of  the 
storm  cases.   The  radially  cumulative  regression  equations  for  the  other 
cases  also  showed  that,  at  the  higher  energy  radii  of  a  storm  (viz, 
R  ^  2°),  larger  inward  (more  negative)  low  cloud  band  crossing  angles  are 
associated  with  higher  windspeed,  while  at  R  >  2°,  smaller  inward  (less 
negative)  cloud  band  crossing  angles  are  associated  with  higher  windspeed. 


The  results  of  the  regression  calculations  for  combined  storm  cases,  by 
radius,  are  shown  in  tables  4  and  5.   Comparison  of  the  two  tables  shows 
the  differences  between  the  algebraic  and  absolute  values  of  O^id,  As 
expected,  the  mean  ahsolute  values  of  ccc-j^  were  larger  than  the  mean 
algebraic  values;  in  either  case,  they  are  poorly  related  to  the  windspeed. 
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Table  4. — Linear  regression  of  observed  windspeed  (vQ^s)  on  algebraic 
values  of  the  cloud  band  crossing  angle  (aci(j)  at  separate  radii  for 
10  combined  cases. 


n 

R  (deg) 

vobs<kt> 

acld  <de8> 

r 

a  (kt) 

b  (kt/deg) 

S  (kt) 

160 

1 

50.0 

-12.84 

-0.17 

47.17 

-0.22 

13.23 

156 

2 

34.6 

-17.93 

-0.10 

32.94 

-0.09 

10.39 

153 

3 

27.4 

-21.42 

0.04 

27.99 

0.03 

8.55 

137 

4 

23.4 

-23.38 

0.21 

25.56 

0.09 

7.48 

123 

5 

20.0 

-25.53 

0.10 

20.82 

0.03 

6.26 

110 

6 

17.6 

-25.84 

0.09 

18.17 

0.02 

5.86 

86 

7 

15.9 

-27.73 

0.16 

16.67 

0.03 

5.32 

64 

8 

14.7 

-30.74 

0.10 

15.31 

0.02 

5.30 

Table  5. — Linear  regression  of  observed  windspeed  (v0^s)  on  absolute 
values  of  the  cloud  band  crossing  angle  (|otc^(j|)  at  separate  radii 
for  10  combined  cases. 


n 

R  (deg) 

Wkt> 

lacldl 

(deg)   r 

a  (kt) 

b  (kt/deg) 

S  (kt 

160 

1 

50.0 

13.53 

0.18 

46.56 

0.25 

13.20 

156 

2 

34.6 

18.25 

0.11 

32.59 

0.11 

10.37 

153 

3 

27.4 

21.52 

-0.03 

27.92 

-0.02 

8.55 

137 

4 

23.4 

25.42 

-0.22 

26.38 

-0.12 

7.47 

123 

5 

20.0 

28.05 

-0.25 

22.65 

-0.09 

6.10 

110 

6 

17.6 

29.66 

-0.14 

18.97 

-0.05 

5.83 

86 

7 

15.9 

34.42 

-0.16 

17.08 

-0.03 

5.32 

64 

8 

14.7 

38.95 

-0.09 

15.74 

-0.03 

5.31 

The  radially  cumulative  regression  results  for  the  10  combined  storm 
cases  are  shown  in  tables  6  and  7. 

When  the  data  from  all  10  cases  were  combined  using  both  algebraic  and 
absolute  values  of  acid  at  separate  radii  (tables  4  and  5)  and  cumulative 
radii  (tables  6  and  7),  the  correlation  changed  sign  at  R  >  3°,  as  in  the 
other  cases.   Although  the  correlation  was  quite  low,  the  relation  shown 
for  the  inner  radii,  (tables  2  to  7)  is  opposite  the  relation  implied  by 
eq.  (10).   The  results  for  the  inner  radii  tend  to  support  the  relation 
derived  by  Malkus  and  Riehl  (1959).   In  that  study,  a  dynamic  model  of  a 
steady-state  mature  hurricane  was  evolved  in  which  the  windspeed  at  dif- 
ferent radii  of  a  storm  was  directly  proportional  to  the  inflow  angle  of 
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Table  6. — Radially  cumulative  linear  regression  of  observed  windspeed 
(v  v  )  on  algebraic  values  of  the  cloud  band  crossing  angle  (01^^) 
for  10  combined  cases. 


n 

2R   (deg) 

obs 

acld    (deg) 

r 
-0.17 

a    (kt) 

b   (kt/deg) 

S    (kt) 

160 

1-1 

50.0 

-12.84 

47.17 

-0.22 

13.23 

316 

1-2 

42.4 

-15.35 

0.01 

42.56 

0.01 

14.29 

469 

1-3 

37.5 

-17.33 

0.13 

40.32 

0.16 

14.39 

606 

1-4 

34.3 

-18.70 

0.20 

38.33 

0.22 

14.25 

729 

1-5 

31.9 

-19.85 

0.22 

36.13 

0.21 

14.17 

839 

1-6 

30.0 

-20.64 

0.22 

34.04 

0.20 

14.17 

925 

1-7 

28.7 

-21.30 

0.22 

32.34 

0.17 

14.17 

989 

1-8 

27.8 

-21.91 

0.22 

31.42 

0.17 

14.15 

Table  7. — Radially  cumulative  linear  regression  of  observed  windspeed 
(v0bs)  on  absolute  values  of  the  cloud  band  crossing  angle  (la^dl) 
for  10  combined  cases. 


n 

R  (deg) 

Wkt> 

Kld|    (deg) 

r 

0.18 

a   (kt) 

b    (kt/deg) 

S    (kt) 

160 

1-1 

50.0 

13.53 

46.56 

0.25 

13.20 

316 

1-2 

42.4 

15.86 

0.00 

42.37 

0.00 

14.29 

469 

1-3 

37.5 

17.71 

-0.12 

40.32 

-0.16 

14.40 

606 

1-4 

34.3 

19.45 

-0.22 

39.44 

-0.26 

14.17 

729 

1-5 

31.9 

20.90 

-0.28 

38.18 

-0.30 

13.95 

839 

1-6 

30.0 

22.05 

-0.30 

36.68 

-0.30 

13.85 

925 

1-7 

28.7 

23.20 

-0.32 

35.32 

-0.29 

13.77 

989 

1-8 

27.8 

24.22 

-0.34 

34.96 

-0.30 

13.65 

air  trajectories,  while  the  relation  shown  in  eq.  (10)  is  based  on  inflow 
angles  of  streamlines. 

The  slight  improvement  in  the  correlation,  when  calculated  for  the  entire 
area,  appears  to  be  largely  due  to  the  relation  between  tangential  windspeed 
and  radial  distance  discussed  in  chapter  I.   (See  eq.  1.)  The  dependence  of 
cloud  band  crossing  angles  on  radial  distance  from  the  storm  center  will  be 
discussed  in  section  III.E. 
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2.   Relative  Windspeed  vs.  Cloud  Band  Crossing  Angle 

Relative  windspeed  (vrei)  is  defined  as  the  magnitude  of  the  vector 
obtained  by  subtracting  the  storm  motion  vector  from  the  observed  wind 
vector.  Coefficients  of  the  linear  regression  of  vre^  on  otc^j  are  shown  in 
tables  8,  9,  10,  and  11.   Comparing  these  tables  with  tables  4,  5,  6,  and 
7,  it  is  seen  that  the  relative  windspeed  was  only  slightly  better  corre- 
lated with  the  cloud  band  crossing  angle  than  was  the  observed  wind. 

Table  8. — Linear  regression  of  relative  windspeed  (vrei)  on  algebraic 
values  of  the  cloud  band  crossing  angle  (acid)  at  separate  radii  for  10 
combined  cases. 


n 

R  (deg) 

vrel(kt) 

acld(deg) 

r 
-0.08 

a  (kt) 

b  (kt/deg) 

S  (kt) 

160 

1 

49.9 

-12.84 

48.64 

-0.10 

12.36 

156 

2 

34.5 

-17.93 

-0.04 

34.01 

-0.03 

8.78 

153 

.3 

27.3 

-21.42 

0.14 

29.41 

0.10 

8.19 

137 

4 

22.7 

-23.38 

0.16 

24.45 

0.08 

7.76 

123 

5 

18.6 

-25.53 

0.15 

20.02 

0.06 

7.68 

110 

6 

16.0 

-25.84 

0.14 

17.22 

0.05 

7.65 

86 

7 

13.8 

-27.73 

0.18 

14.93 

0.04 

7.44 

64 

8 

12.8 

-30.74 

0,19 

14.37 

0.05 

7.84 

Table  9. — Linear  regression  of  relative  windspeed  (vrei)  on  absolute 
values  of  the  cloud  band  crossing  angle  C l^cld I ^  at  seParate  radii 
for  10  combined  cases. 


n 

R  (deg) 

vrel(kt) 

lacldl <de8> 

r 
0.10 

a  (kt) 

b  (kt/deg) 

S  (kt) 

160 

1 

49.9 

13.53 

48.13 

0.13 

12.34 

156 

2 

34.5 

18.25 

0.04 

33.85 

0.04 

8.78 

153 

3 

27.3 

21.52 

-0.14 

29.49 

-0.10 

8.19 

137 

4 

22.7 

25.42 

-0.27 

26.57 

-0.15 

7.57 

123 

5 

18.6 

28.05 

-0.29 

22.44 

-0.14 

7.42 

110 

6 

16.0 

29.66 

-0.19 

18.41 

-0.08 

7.58 

86 

7 

13.8 

34.42 

-0.06 

14.42 

-0.02 

7.55 

64 

8 

12.8 

38.95 

-0.09 

14.31 

-0.04 

7.96 

As  in  the  observed  winds  case,  the  correlation  coefficients  are  small, 
but  the  change  in  sign  at  radius  3°  is  again  present.   The  regression 
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results  for  separate  radii,  shown  in  tables  8  and  9,  reveal  that  the 
absolute  values  of  acid  were  only  slightly  better  correlated  with  the 
relative  wind  speed  than  were  the  algebraic  values.  This  was  also  true 
for  the  observed  wind.  Aside  from  some  possible  significance  of  the 
change  in  sign  between  2°  and  3°  radii,  it  appears  that  small  inward 
crossing  angles  of  the  cloud  bands  tend  to  accompany  higher  windspeed 
while  large  crossing  angles  tend  to  accompany  low  windspeed.  This  ten- 
dency was  expected,  but  much  higher  correlations  were  anticipated. 

Tables  10  and  11  show  the  results  of  the  radially  cumulative  linear 
regression  between  vrei  and  0t  i  .  for  the  combined  storm  cases.  Table  10 
shows  the  results  obtained  by  using  the  algebraic  values  and  table  11  the 
absolute  values  of  aci<j. 


Table  10. — Radially  cumulative  linear  regression  of  relative  windspeed 
(yrei)  on  algebraic  values  of  the  cloud  band  crossing  angle  (aci(j)  f°r 
10  combined  cases. 


n 

ZR  (deg) 

vrel(kt) 

acld(deg) 

r 
-0.08 

a  (kt) 

b  (kt/deg) 

S  (kt) 

160 

1-1 

49.9 

-12.84 

48.64 

-0.10 

12.36 

316 

1-2 

42.3 

-15.35 

0.08 

43.82 

0.10 

13.19 

469 

1-3 

37.4 

-17.33 

0.20 

41.56 

0.24 

13.49 

606 

1-4 

34.1 

-18.70 

0.24 

38.89 

0.26 

13.67 

729 

1-5 

31.5 

-19.85 

0.26 

36.52 

0.26 

13.95 

839 

1-6 

29.4 

-20.64 

0.26 

34.27 

0.23 

14.21 

925 

1-7 

28.0 

-21.30 

0.26 

32.35 

0.21 

14.43 

989 

1-8 

27.0 

-21.91 

0.26 

31.39 

0.20 

14.51 

Table  11. — Radially  cumulative  linear  regression  of  relative  windspeed 
(vre^)  on  absolute  values  of  the  cloud  band  crossing  angle  (l^ridl)  for 
10  combined  cases. 


n 

ZR  (deg) 

vrel(kt) 

Kldl(deg) 

r 

a  (kt) 

b  (kt/deg) 

S  (kt) 

160 

1-1 

49.9 

13.53 

0.10 

48.13 

0.13 

12.34 

316 

1-2 

42.3 

15.86 

-0.07 

43.74 

-0.09 

13.20 

469 

1-3 

37.4 

17.71 

-0.19 

41.68 

-0.24 

13.52 

606 

1-4 

34.1 

19.45 

-0.29 

40.52 

-0.33 

13.50 

729 

1-5 

31.5 

20.90 

-0.34 

39.13 

-0.37 

13.60 

839 

1-6 

29.4 

22.05 

-0.36 

37.47 

-0.36 

13.74 

925 

1-7 

28,0 

23.20 

-0.36 

35.72 

-0.33 

13.92 

989 

1-8 

27.0 

24.22 

-0.38 

35.29 

-0.34 

13.90 
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The  correlation  coefficients  for  the  entire  data  sets  were  0,26  for  the 
algebraic  values  and  -0,38  for  the  absolute  values.  The  variance  of 
relative  windspeed  explained  was  6.9%  for  the  algebraic  values  of  a  ^.d 
and  14.6%  for  the  absolute  values.   Comparable  percentages  for  the 
observed  windspeed  were  5.0%  and  11.6%,  respectively.   As  in  the  observed 
wind  case,  the  increase  in  the  correlation  as  more  and  more  radii  were 
included  indicates  that  a  -^  is  related  to  vre^  because  aci(j  itself  is 
related  to  the  radial  distance  (R)  from  the  storm  center, 

B.   Linear  Regression  of  Windspeed  on  Infrared  Temperature  (TBB) 

1,   Observed  Windspeed  vs.  Infrared  Temperature 

The  procedure  used  to  test  for  relation  between  the  windspeed  and  the 
cloud  band  crossing  angle  was  repeated,  except  that  infrared  temperature 
was  used  in  place  of  the  cloud  band  crossing  angle.  Table  12  shows  the 
linear  regression  of  v0bs  on  TBB  at  separate  radii  for  the  Sept.  1 
hurricane  Carmen  case. 


Table  12. — Linear  regression  of  observed  windspeed  (v0bs)  on  infrared 
temperature  (TBB)  at  separate  radii  for  hurricane  Carmen,  1600  GMT, 
Sept.  1,  1974. 


n 

R  (deg) 

vobs<kt) 

tbb(°k> 

213.6 

r 
0.11 

a  (kt) 
0.80 

b  (kt/°K) 

S  (kt) 

16 

1 

43.9 

0.20 

11.80 

16 

2 

25.8 

223.1 

-0.35 

63.60 

-0.17 

8.38 

16 

3 

20.4 

245.8 

-0.46 

48.74 

-0.12 

5.90 

16 

4 

16.9 

261.4 

-0.27 

31.90 

-0.06 

5.96 

16 

5 

15.3 

284.4 

-0.52 

70.94 

-0.20 

4.58 

16 

6 

13.6 

286.4 

0.59 

-66.24 

0.28 

4.20 

16 

7 

12.3 

284.3 

0.58 

-41.01 

0.19 

4.06 

16 

8 

11.4 

278.4 

0.07 

7.81 

0.01 

3.91 

For  this  case,  the  TBB  increased  to  a  maximum  at  the  6®  radius  and  then 
slowly  decreased.   The  v0bs  decreased  monotonically  with  increasing  radius, 
The  correlations  at  separate  radii,  in  this  and  the  remaining  cases,  was 
small  and  negative. 

The  radially  cumulative  regression  results  for  the  Sept.  1,  1974 
hurricane  Carmen  case  are  shown  in  table  13. 
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Table  13. — Radially  cumulative  linear  regression  of  observed  wlndspeed 
(v0|jS)  on  infrared  temperature  (TBB)  for  hurricane  Carmen,  1600  GMT, 
Se,pt.  1,  1974. 


n 

ER  (deg) 

vobs(kt) 

TBB(°K) 
213.6 

r 
0.11 

a  (kt) 

b  (kt/°K) 

S  (kt) 

16 

1-1 

43.9 

0.80 

0.20 

11.80 

32 

1-2 

34.8 

218.4 

-0.33 

103.64 

-0.32 

13.07 

48 

1-3 

30.0 

227.5 

-0.50 

97.20 

-0.30 

11.92 

64 

1-4 

26.7 

236.0 

-0.55 

88.19 

-0.26 

11.31 

80 

1-5 

24.4 

245.7 

-0.62 

85.69 

-0.25 

10.33 

96 

1-6 

22.6 

252.4 

-0.64 

84.15 

-0.24 

9.87 

112 

1-7 

21.2 

257,0 

-0.64 

82.83 

-0.24 

9.63 

128 

1-8 

20.0 

259.7 

-0.64 

81.56 

-0.24 

9.43 

The  correlation  coefficients  increased  as  additional  radii  were  included 
in  the  data  set,  as  was  the  case  for  the  cloud  band  crossing  angles. 

As  noted  earlier,  infrared  temperatures  were  available  for  only  8  of  the 
10  storm  cases.  Table  14  shows  the  regression  of  observed  wlndspeed  (v 0hs) 
on  infrared  temperature  (TBB)  at  separate  radii  for  the  8  combined  storm 
cases. 


Table  14. — Linear  regression  of  observed  wlndspeed  (v0bs)  on  infrared 
temperature  (TBB)  at  separate  radii  for  8  combined  cases. 


n 

R  (deg) 

v0bs(kt> 

TBB(°K) 

r 
-0.23 

a  (kt) 

b  (kt/°K) 

S  (kt) 

128 

1 

50.4 

216.5 

87.58 

-0.17 

11.40 

128 

2 

35.7 

232.6 

-0.25 

62.94 

-0.12 

9.71 

128 

3 

28.3 

247.8 

-0.18 

43.52 

-0.06 

8.39 

126 

4 

23.6 

253.5 

-0.14 

32.80 

-0.04 

7.48 

124 

5 

19.5 

265.0 

-0.18 

31.04 

-0.04 

6.79 

124 

6 

16.3 

272.9 

-0.09 

22.53 

-0.02 

6.00 

123 

7 

13.7 

276.2 

0.02 

12.20 

0.01 

5.32 

122 

8 

12.6 

280.2 

-0.01 

13.71 

-0.00 

5.06 

Small,  negative  values  of  r  show  that  infrared  temperature  is  poorly 
related  to  the  wlndspeed. 

Table  15  shows  the  radially  cumulative  linear  regression  of  observed 

wlndspeed  (v0fos)  on  infrared  temperature  (TBB)  for  8  combined  storm  cases, 
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Table  15, — Radially  cumulative  linear  regression  of  observed  windspeed 
(v0bs)  on  infrared  temperature  (Tug)  for  8  combined  cases. 


n 

£R  (deg) 

vobs<kt> 

TBB('K) 

r 
-0.23 

a  (kt) 

b  (kt/°K) 

S  (kt) 

128 

1-1 

50.4 

216.5 

87.58 

-0.17 

11.40 

256 

1-2 

43.0 

224.5 

-0.40 

100,96 

-0.26 

12.06 

384 

1-3 

38.1 

232.3 

-0.48 

99.97 

-0.27 

12.05 

510 

1-4 

34.5 

237.5 

-0.49 

93.50 

-0.25 

12.19 

634 

1-5 

31.6 

242.9 

-0.53 

93.08 

-0.25 

12.07 

758 

1-6 

29.1 

247.8 

-0.56 

94.25 

-0.26 

11.87 

881 

1-7 

26.9 

251.8 

-0.59 

95.32 

-0.27 

11.73 

1003 

1-8 

25.2 

255.2 

-0.61 

96.87 

-0.28 

11.47 

For  the  entire  data  set,  37.2%  of  the  observed  windspeed  variance  could 
be  explained  by  the  infrared  temperature;  however,  the  low  correlation 
coefficients  calculated  at  each  separate  radius  (table  14)  indicate  that, 
as  in  the  case  of  the  cloud  band  crossing  angle,  the  relation  clearly 
depends  on  the  radial  distance  from  the  storm  center. 

2.   Relative  Windspeed  vs.  Infrared  Temperature 


To  determine  whether  relative  windspeed  was  better  related  to  infrared 
temperature  than  was  observed  windspeed,  linear  regressions  were  again 
computed  using  relative  windspeed  as  the  dependent  variable.  The  regres- 
sions showed,  as  in  the  observed  wind  case,  that  at  individual  radii  the 
correlation  was  low  and  fluctuated  in  sign.  When  the  data  for  all  radii 
were  combined,  the  correlation  increased  significantly.   For  the  entire 
data  set,  the  correlation  was  -0.595.  Thus,  about  35.4%  of  the  variance 
of  the  relative  windspeed,  as  opposed  to  37.2%  of  the  variance  of  the 
observed  windspeed,  could  be  accounted  for  by  the  infrared  temperature. 
The  1.8%  difference  in  the  explained  variance  is  considered  insignificant. 

As  noted  earlier  in  the  observed  wind  case,  the  relation  of  infrared 
temperature  to  relative  windspeed  also  existed  primarily  because  the 
infrared  temperature  was  related  to  radial  distance  from  the  storm  center. 
The  predominance  of  negative  correlation  coefficients  at  separate  radii 
indicated  that  the  expected  inverse  relation  between  the  windspeed  and 
infrared  temperature  exists,  but  is  too  weak  to  be  useful.   The  use  of 
infrared  temperature  to  specify  the  areal  windspeed  distribution  will  be 
discussed  further  in  section  III.  E. 
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C,  Linear  Regression  of  Cloud  Band  Crossing  Angle  on 
Wind  Field  Crossing  Angle 

1.   Cloud  Band  Crossing  Angle  vs.  Observed  Wind  Crossing  Angle 

The  cloud  band  crossing  angle  was  selected  as  one  of  the  satellite  data 
parameters  because  it  was  assumed  that  this  angle  would  be  well  correlated 
with  the  wind  crossing  angle  (ct0t>s) ,  which  in  turn  was  shown  In  eq,  (10) 
to  be  related  to  the  tangential  component  of  the  windspeed. 

How  well  the  crossing  angles  are  correlated  indicates  how  close  the  wind 
direction  is  to  the  orientation  of  low-level  cloud  bands.   Table  16  shows 
results  of  the  linear  regression  of  algebraic  values  of  otc^(j  on  algebraic 
values  of  the  crossing  angle  of  the  wind  direction  (ot  ^s)  at  separate  radii 
for  the  10  combined  storm  cases. 


Table  16. — Linear  regression  of  algebraic  values  of  the  cloud  band  crossing 
angle  (cx^d)  on  algebraic  values  of  the  wind  field  crossing  angle  (ct0^s) 
at  separate  radii  for  10  combined  cases. 


n 

R  (deg) 

acld(deg) 

aobs(de8> 

r 
0.00 

a  (deg) 
-12.84 

b 
0.00 

S  (deg) 

160 

1 

-12.84 

-  6.46 

10.68 

156 

2 

-17.93 

-  8.57 

0.19 

-16.90 

0.12 

10.82 

153 

3 

-21.42 

-11.79 

0.38 

-18.94 

0.21 

10.79 

137 

4 

-23.38 

-18.42 

0.43 

-17,44 

0.32 

15.47 

123 

5 

-25.53 

-25.85 

0.42 

-18.16 

0.29 

18.33 

110 

6 

-25.84 

-31.07 

0.38 

-18.29 

0.24 

21.72 

86 

7 

-27.73 

-31.42 

0.32 

-20.59 

0.23 

30.38 

64 

8 

-30.74 

-29.24 

0.48 

-21.78 

0.31 

26.78 

Table  17  shows  the  results  obtained  by  using  the  absolute  values  of  the 
crossing  angles  of  both  the  cloud  bands  and  the  wind  direction.   This 
table  shows  the  magnitudes  of  the  crossing  angles  independent  of  their 
inward  (-)  or  outward  (+)  orientation. 
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Table  17. — Linear  regression  of  absolute  values  of  the  cloud  band  crossing 
angle  (|otc^^|)  on  absolute  values  of  the  wind  field  crossing  angle 
(|a0bs|)  at  separate  radii  for  10  combined  cases. 


p 

R  (deg) 

hcldlWeg) 

laobsl(de8> 

r 
0.08 

a   (deg) 

b 
0.70 

S   (deg) 

160 

I 

13.53 

13.47 

12.59 

10.65 

156 

2 

18.25 

16.01 

0.16 

15.68 

0.16 

10.88 

153 

3 

21.52 

20.04 

0.34 

15.38 

0.31 

10.94 

137 

4 

25.42 

25.20 

0.30 

16.84 

0.34 

16.35 

123 

5 

28,05 

31.97 

0.26 

20.66 

0.23 

19.47 

110 

6 

29.66 

39.94 

0.24 

21.24 

0.21 

22.80 

86 

7 

34.42 

46.99 

0.12 

27.97 

0.14 

31.83 

64 

8 

38.95 

49.91 

0.28 

23.01 

0.32 

29.41 

Comparison  of  the  algebraic  means  with  the  absolute  means  at  each  radius 
in  tables  16  and  17  reveals  the  extent  to  which  the  outward-oriented 
crossing  angles  (+  values)  reduced  the  average  inward  algebraic  crossing 
angles  (-  values)  at  each  radius.   Table  17  shows  that  the  mean  absolute 
values  of  the  crossing  angle  of  the  observed  wind  (|0iobs|)  became  larger 
than  the  mean  absolute  values  of  the  cloud  band  crossing  angle  C ! OtcldT^ 
beyond  the  4°  radius.   This  result  was  anticipated  since  the  upper  limits 
of  the  cloud  band  crossing  angles  were  +  90°,  while  the  upper  limits  for 
the  observed  wind  crossing  angles  were  +  180°. 

It  was  expected  that  the  correlation  coefficients  would  be  higher  at  all 
radii  than  those  shown  in  table  16,  and  that  the  correlations  near  the 
storm  center  would  be  higher  than  those  found  at  the  outer  radii.   Neither 
expectation  was  fulfilled.   Figure  9  illustrates  the  ranges  of,  and  the 
lack  of  correlation  between,  the  crossing  angle  of  the  wind  direction  (a0t>s) 
and  the  crossing  angle  of  the  low  cloud  bands  (acid)  for  the  1°  radius 
for  the  10  combined  storm  cases.   The  lack  of  correlation  can  also  be 
seen  by  simple  inspection  of  figure  9. 
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Figure  9. — Wind  field  crossing  angle  (ot0bs)  vs.  cloud  band  crossing  angle 
0*cld)  for  the  1°  radius  for  10  combined  storm  cases.   (N  -  160). 
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To  what  extent  the  results  presented  In  table  16  are  representative  is 
largely  dependent  upon  the  reliability  of  the  wind  direction.  Even  though 
a  field  analysis  method  was  used  to  extract  the  wind  directions  and  speeds, 
the  areal  gaps  in  the  observations  may  have  been  too  large,  or  the  +  8-hour 
compositing  time  interval  may  have  been  too  long  to  yield  an  accurate  wind 
field.  Some  of  the  poor  correlation  may  also  have  been  caused  by  differ- 
ences in  altitude  of  the  winds  and  the  cloud  bands,  particularly  in  areas 
with  large  vertical  shear  in  the  wind  direction. 

Because  of  the  uncertainties  in  the  wind  fields,  and  the  limited  number 
of  storm  cases,  only  tentative  conclusions  can  be  drawn. 

2,   Cloud  Band  Crossing  Angle  vs.  Relative  Wind  Crossing  Angle 

Regressions  similar  to  those  shown  in  tables  16  and  17  were  computed 
using  the  relative  wind  crossing  angle  (arei) .  The  correlations  (not 
shown)  of  the  cloud  band  crossing  angle  with  the  relative  wind  crossing 
angle  were  even  lower  than  for  the  observed  wind  crossing  angle.  The 
highest  correlation  coefficient  was  0.24  at  the  3°  radius,  followed  by  0.14 
at  5%  and  0.13  at  4°. 

D.  Linear  Regression  of  Windspeed  on  Wind  Crossing  Angle 

In  section  III. A,  it  was  shown  that  v0|jS  and  clci^   at  individual  radii 
were  poorly  correlated.  When  all  radii  of  the  storms  were  combined,  thereby 
introducing  the  relation  of  vODS  to  R,  the  correlation  was  still  only  about 
0,4.  Further,  since  in  section  III.C,  it  was  shown  that  the  correlation 
between  aci^  and  a0^,s  at  all  radii  was  below  0.5,  the  question  arose  as  to 
whether  a0bs  itself  was  related  to  the  windspeed.  Accordingly,  the  linear 
regression  of  v0^s  on  a  i  was  examined.  Similar  computations  were  also 
made  for  the  relative  wind  field. 

At  separate  radii  neither  a0^s  nor  ctrei   showed  significant  Improvement 
over  ac^  in  their  correlation  to  the  windspeed.  Cumulatively  for  all 
radii,  the  algebraic  values  of  0to^s  and  otre^  showed  a  correlation  of  only 
about  0.1  to  the  windspeed  while  the  absolute  values  of  a0^,s  showed  a 
correlation  of  about  0.4  with  windspeed.  Comparative  values  for  the  cloud 
band  crossing  angle  were  0.25  for  algebraic  values  and  0.36  for  absolute 
values.  Overall,  the  windspeed  was  slightly  better  correlated  with  aci(j 
than  cx0fos.  This  also  was  unexpected. 

E.  Screening  Regression 

1.   Procedure 

The  radially  cumulative  correlation  coefficients  presented  earlier 
increased  monotonlcally  as  data  from  additional  radii  were  included. 
Apparently,  the  correlation  of  these  two  satellite  parameters  with  windspeed 
is  related  to  the  radial  distance  from  the  storm  center,  as  shown  by 
eq.  (1). 
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To  evaluate  the  radial  relation  and  the  effectiveness  of  the  two  satel- 
lite data  parameters  in  specifying  the  wlndspeed  distribution,  a  Stagewise 
Screening  Regression  Method*  was  applied  to  the  data.  Only  selected 
results  of  this  analysis  will  be  presented. 

Later  in  the  study,  a  more  rigorous  statistical  analysis,  the  Stepwise 
Multiple  Regression  Method,*  was  applied  to  the  same  data.  The  results 
obtained  by  this  method  showed  only  minor  differences  from  those  obtained 
earlier  with  the  Stagewise  Screening  Method,  thereby  corroborating  the 
results  of  the  Stagewise  Screening  Method.  In  both  Screening  Regression 
Methods,  any  number  of  variables  can  be  entered  as  predictors.  Further, 
nonlinear  operations  can  be  performed  on  any  of  the  variables  before  the 
regression  process  has  begun. 

In  applying  screening  regression  methods,  the  wlndspeed  was  usually  desig- 
nated as  the  dependent  variable.  However,  separate  runs  with  the  Stagewise 
Screening  Method  were  made  by  using  first  the  radial  and  then  the 
tangential  variation  of  the  wlndspeed  as  the  dependent  variable.  The 
principal  Independent  variables  tested  were:   the  radius  (R) ,  the  cloud 
band  crossing  angle  (a) ,  the  infrared  temperature  (Tjjg) ,  and  the  radial  and 
tangential  variation  of  a  and  Tgg.  The  9  forms  of  the  variables  are  shown 
later  in  column  7  of  table  18. 

2,   Screening  Regression  of  Observed  Wlndspeed  on  Nine  Independent 
Variables 

Because  coincident  values  of  each  of  the  variables  are  required  for  each 
grid  point,  the  sample  was  reduced  significantly.  A  further  reduction  of 
the  sample  occurred  because  spatial  differences  of  the  variables,  which 
require  simultaneous  data  at  two  points,  were  also  used  as  variables. 

For  the  variables  screened  by  both  methods,  as  shown  in  table  18,  the 
absences  of  1  or  more  of  the  9  independent  variables  reduced  the  sample 
from  a  maximum  possible  1120  to  601.  Because  2  additional  samples  had  been 
retained  in  the  Stepwise  method,  the  sample  size  then  was  603. 

The  left  side  of  table  18  shows  the  results  of  the  Stagewise  Screening 
Method.  Only  the  first  10  passes  through  the  data  are  shown  since  the 
variance  explained  by  additional  passes  would  be  negligible.  Under  the 
Stepwise  Multiple  Regression,  the  step  numbers  show  the  order  of  selection 
followed  by  the  9  forms  of  the  independent  variables  screened. 


*See  appendix  for  a  description  of  screening  regression  methods. 
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Table  18. — Results  of  both  the  Stagewise  Screening  and  the  Stepwise 

Multiple  Regressions  of  the  observed  windspeed  (dependent  variable)  on 
9  independent  variables  for  all  10  storm  cases. 


Stagewise  Screening  Regression 

N  -  601 

Indiv.   Var. 
Pass  Cor.    Expl. 

no.   Variable  Coef.     %    A% 


Stepwise  Multiple  Regression 

N  =  603 

Cum.   Var. 
Step  Cor.   Expl. 

no.  Variable  Coef.   %     A% 


1* 

R-0.5 

.7932 

62.9 

0.4 

1 

R-0.5 

.7912 

62.61 

0.59 

2 

Aa/AR 

.0988 

63.3 

0.4 

2 

TBB 

.7950 

63.20 

0.47 

3+ 

TBB 

-.1080 

63.7 

0.3 

3 

Aa/AR 

.7979 

63.67 

0.36 

4 

ATBB/AR 

-.0835 

64.0 

0.2 

4 

ATBB/A0 

.8002 

64.03 

0.27 

5 

Aa/AG 

.0774 

64.2 

0.1 

5 

Aa/AG 

.8019 

64.30 

0.23 

6** 

R-0.5 

-.0418 

64.3 

0.1 

6 

a 

.8033 

64.53 

0.23 

7++ 

TBB 

-.0646 

64.4 

0.1 

7 

ATBB/AR 

.8048 

64.76 

0.06 

8 

ATBB/AG 

-.0582 

64.5 

0.1 

8 

R 

.8051 

64.82 

0.05 

9 

a 

.0430 

64.6 

0.1 

9 

|a| 

.8054 

64.87 

10*** 

R-0.5 

-.0554 

For  the  Stagewise  Screening  method,  where  alternating  selection  of  the 
same  variable  is  permitted,  R~0«5  was  selected  on  the  first,  sixth,  and 
_enth  passes,  TBB  was  selected  on  the  third  and  seventh  passes,  while  the 
absolute  value  of  the  cloud  band  crossing  angle  (|a|)  and  the  linear  radius 
(R)  were  never  selected. 

Both  the  Stepwise  and  Stagewise  methods  show  correlation  coefficients  of 
.0,79,  i.e.,  63%  of  the  variance  of  the  windspeed  could  be  explained  by  the 
radius.   The  first  pass/step  of  the  two  methods  is  mathematically  identi- 
cal, but  in  subsequent  passes  the  procedures  differ  (see  appendix). 
However,  both  methods  selected  either  the  infrared  temperature  (Tbb)  or  tlie 
radial  variation  of  the  cloud  band  crossing  angle  (Aa/AR)  as  the  second 
and  third  variables  in  order  of  their  relation  to  the  windspeed.  Note  that 
the  correlation  coefficients  are  given  individually  in  the  Stagewise  but 
cumulatively  in  the  Stepwise  Method. 
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By  either  method,  only  about  2%  additional  variance  of  the  windspeed  can 
be  accounted  for  by  TBB  and  o;,  above  the  63%  already  accounted  for  by  R~0.5, 
The  regression  equations  derived  for  each  method  are  shown  in  table  18a, 


Table  18a. — Screening  Regression  Equations  derived  for  the  observed 

windspeed  (v0^s)  and  the  9  independent  variables  shown  in  table  18,  using 
Stagewise  and  Stepwise  methods. 

Stagewlse  equation:  vobs  =  12.27  +  51.01[R-°«5]  +  .09[Act/AR]  -  .05|TBB] 
-  .03[ATBB/ARJ  +  .O4[Aa/A0]  -  .O31ATBB/A0]  +  .02 [a] 

Stepwise  equation:  v0t>s  =  27.12  +  42.85[R~0'5]  -  .08[TBB]  +  ,09[Aa/AR] 
-  .04[ATBB/AR]  +  .O8[Aa/A0J  +  .08[aJ  -  .04[ATBB/AR]  -  .60[R]  +  .041 |a|] 


These  equations  represent  the  best  fit  between  the  windspeed  and  the 
nine  independent  variables  screened.   The  coefficients  indicate  the 
relative  importance  of  each  variable  only  when  the  proper  units  are 
applied.   The  radial  increments  (AR)  and  the  tangential  increments  (AG) 
were  treated  as  unit  values  since  they  were  always  taken  as  1°  latitude 
in  the  radial  direction  and  22.5°  (ir/8  radians)  in  the  tangential  direction. 

3.   Individual  Correlations  between  Windspeed  and  Nine  Independent  Variables 


Because  only  the  variable  with  the  highest  correlation  was  shown  in  the 
first  pass  of  the  screening  methods,  the  individual  correlations  of  the 
remaining  variables  are  not  shown  in  table  18.   Table  19  shows  the  indi- 
vidual relations  between  the  windspeed  and  each  of  the  nine  independent 
variables.   It  is  instructive  to  evaluate  these  correlation  coefficients 
(keeping  in  mind  that  they  are  not  statistically  independent)  to  see  which 
variables,  or  forms  of  the  variables,  are  most  closely  related  to  the 
windspeed . 

Table  19. — Correlation  (r)  between  the  windspeed  (v)  and  the  9  independent 
variables  shown  in  table  18.   (N  =  601). 


No. 

Variable 

r 

No. 

Variable 

r 
-.3523 

No. 
7 

Variable 

r 

1 

R-0.5 

.7932 

4 

|  ex  | 

TBB/A0 

-.0594 

2 

R 

-.7473 

5 

a 

.2461 

8 

Aa/A0 

.0203 

3 

T 
XBB 

-.5383 

6 

ATBB/AR 

.0934 

9 

Aa/AR 

-.0071 
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In  table  19,  the  correlation  coefficients  of  the  nine  independent  vari- 
ables are  listed  in  decreasing  order  of  importance.   The  subscripts  "eld" 
and  "obs"  have  been  omitted  for  convenience.   Both  forms  of  the  radius  (R) 
were  reasonably  well  correlated  with  the  windspeed.  However,  the  radial 
and  tangential  variations  of  TgB  and  a^  rank  at  the  bottom  with  corre- 
lation coefficients  less  than  0.1. 


Because  eq.  (10)  suggested  the  cloud  band  crossing  angle  as  one  of  the 
satellite  data  parameters,  a  combined  form  of  the  variable,  [Cos  acld/R]0,5 
was  selected  for  screening.   Table  20  shows  the  first  pass  of  the  Stagewise 
Screening  Regression  run  that  included  this  variable. 


Table  20. — Correlation  (r)  between  the  windspeed  (v)  and  the  9  independent 
variables  listed  below.   (N  =  717) . 


No. 

Variable      r 

No. 

Variable 

r 

No. 

Variable 

r 

1 

R-0-5     .8021 

4 

TBB 

-.5553 

7 

(Cos  a)0'5 

.3299 

2 

,   -0.5 
(Cosot/R)     .7929 

5 

Cos  a 

.3428 

8 

a 

.2325 

3 

R      -.7539 

6 

|a| 

-.3365 

9 

0 

-.2243 

In  table  20,  the  variables  are  listed  in  decreasing  order  of  their 
relation  to  windspeed.   Note  that  the  variable  R~0.5  is  again  best  corre- 
lated to  the  windspeed,  with  the  combined  form  a  close  second.  The  combined 
form  of  the  variable  showed  no  improvement  over  R~0'5  itself. 

Because  the  storm  cases  in  this  study  were  of  different  intensities,  the 
differences  in  windspeed  might  have  affected  the  relations.   In  an  attempt 
to  isolate  such  a  bias,  as  well  as  to  test  whether  the  variation  of  the 
windspeed  in  the  radial  or  tangential  directions  was  related  to  the  other 
variables,  the  radial  (Av/AR)  and  tangential  (Av/AO)  variations  of  the 
windspeed  were  considered  separately  as  the  dependent  variable. 

4.   Individual  Correlations  between  Spatial  Variations  of  the  Windspeed 
and  Nine  Independent  Variables 

Tables  21  and  22  show  the  correlation  between  spatial  variations  of  the 
windspeed  and  the  9  variables  shown  in  table  18.  The  coefficients  in 
table  21  show  the  importance  of  the  independent  variables  in  decreasing 
order.  R~0«5  was  best  related  to  Av/AR;  while  again,  the  spatial  varia- 
tions rank  at  the  bottom  of  the  list. 
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Table  21. — Correlation  (r)  between  the  radial  variation  of  the  windspeed 
(Av/AR)  and  the  9  independent  variables  shown  in  table  18.   (N  =  601). 


No. 

Variable 

r 
-.6495 

No. 
4 

Variable 

r 

No. 
7 

Variable 

r 

1 

R-0.5 

i«l 

.2448 

Aa/AR 

.0628 

2 

R 

.5788 

5 

a 

-.1563 

8 

ATBB/A0 

.0436 

3 

TBB 

.4232 

6 

ATBB/AR 

-.0892 

9 

Aa/A0 

-.0340 

The  results  of  table  21  show  only  the  first  pass  of  the  Stagewise 
Screening  Regression.   However,  when  10  passes  were  made  through  the  data 
with  (Av/AR)  as  the  dependent  variable,  the  variance  of  the  radial  varia- 
tion of  the  windspeed  is  reduced  by  only  43.2%,  compared  with  64.7%  shown 
in  table  18.   Based  on  the  results  here,  the  radial  variations  of  the 
windspeed  show  a  lower  correlation  on  six  of  the  nine  variables  when 
compared  to  those  obtained  for  the  windspeed  itself. 

Table  22  gives  the  results  of  the  first  pass  of  the  Stagewise  Screening 
Regression  in  which  the  tangential  variation  of  the  windspeed  was  used  as 
the  dependent  variable. 


Table  22. — Correlation  (r)  between  the  tangential  variation  of  the  wind- 
speed  (Av/A0)  and  the  9  independent  variables  shown  in  table  18. 
(N  =  601). 


No. 

Variable 

r 

No. 
4 

Variable 

r 

No. 
7 

Variable 

r 

1 

a 

-.1300 

R-0.5 

-.0435 

ATBB/A0 

.0185 

2 

Aa/AR 

-.1061 

5 

R 

.0358 

8 

Aa/A0 

.0114 

3 

l«l 

.1021 

6 

TBB 

.0260 

9 

ATBB/AR 

-.0050 
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In  table  22,  all  9  independent  variables  show  a  poor  relation  to  the 
tangential  variation  of  windspeed;  however,  R~0«5  has  now  been  replaced  by 
cloud  band  crossing  angle  (a) . 

After  10  screening  passes,  3.7%  reduction  in  variance  of  Av/A0  was 
attained  (result  not  shown) .   This  small  amount  of  explained  variance  in 
the  windspeed  is  far  below  what  might  be  considered  operationally  useful; 
however,  3.3%,  or  nearly  all  the  explained  variance,  was  accounted  for  by 
the  variables  a  and  Aa/AR.  The  negative  signs  of  the  correlation  coeffi- 
cients for  these  variables  were  expected  from  dynamical  considerations. 
Since  a  was  designated  negative  when  inward  (implying  confluence) ,  the 
negative  correlation  coefficient  for  a  indicates  that,  in  the  downstream 
counterclockwise  direction  interval  AG,  Av  would  be  higher  if  the  crossing 
angles  were  more  inward.   Conversely,  if  the  crossing  angles  were  outward 
(Implying  dif luence) ,  the  windspeed  would  decrease  downstream.   Similar 
changes  of  Av/AG  are  indicated  by  evaluating  the  confluence/difluence 
implied  by  the  sign  of  Aa/AR. 

Finally,  it  was  found  that  the  radial  variation  of  v0bs  was  not  as 
closely  related  to  the  independent  variables  as  was  the  windspeed  itself. 
Thus,  any  bias  that  might  have  been  introduced  by  combining  storms  of 
different  intensities  was  not  detected. 

From  the  results  of  the  Screening  Regression,  about  63%  of  the  variance 
of  the  windspeed  can  be  explained  by  R~0*->.  An  additional  2%  to  3%  of  the 
variance  of  the  windspeed  in  the  radial  direction  can  be  accounted  for  by 
the  various  forms  of  Tgg  and  clc^,    and  as  much  as  3%  of  the  variance  in 
the  tangential  direction  can  be  accounted  for  by  ac±d   and  its  radial 
variation. 

IV.   SUMMARY 

Ten  storm  cases  were  studied  to  determine  whether  the  areal  windspeed 
distribution  could  be  estimated  with  the  use  of  two  satellite  data  param- 
eters, the  cloud  band  crossing  angle  (otj,^)  and  the  infrared  temperature 
(Tbb) •   The  correlations  of  these  parameters  with  windspeed  over  the  entire 
storm  area  were  approximately  0.4  for  a.c±d   and  -0.6  for  Tgg.  However, 
these  two  parameters  were  found  to  be  related  to  the  radial  distance  from 
the  storm  center,  R;  and  since  it  was  known  from  earlier  studies  that  the 
tangential  component  of  the  windspeed  can  be  approximated  by  R~0«5,  the 
two  parameters  would  be  useful  only  to  the  extent  that  they  are  statis- 
tically independent  of  R. 

When  linear  correlations  of  windspeed  with  the  satellite  data  parameters 
were  calculated  for  each  storm  case  for  individual  radii,  the  correlation 
coefficients  were  low,  unstable,  and  even  changed  sign  from  radius  to 
radius.   Further,  for  all  10  storm  cases  combined,  both  the  algebraic  and 
absolute  values  of  olci^   were  poorly  correlated  with  both  the  observed  (v0bs) 
and  relative  (vre^)  windspeed  at  separate  radii.  Most  of  the  correlation 
coefficients  were  below  0.2.   The  correlations  between  Tgg  and  v0bs  and 
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vrel  a"t  separate  radii  were  also  smaller  than  0.2,  (See  table  14.).   Beyond 
the  2°  radius,  the  signs  of  the  coefficients  were  as  expected  (See  section 
II. A),  i.e.,  the  correlations  between  the  windspeed  and  the  algebraic  value 
of  acid  were  mostly  positive  (e.g.,  table  4).  The  positive  coefficients 
result  because  the  windspeed  decreased  as  add  became  more  inward  (more 
negative).   The  only  exception  occurred  at  the  innermost  radii  (R  =  1°  and  2°) 
where  the  coefficients  were  negative. 

Relative  windspeed  was  only  slightly  better  correlated  with  ac^d  tnan 
was  the  observed  windspeed.   Further,  the  absolute  value  of  <Xci^   was  only 
slightly  better  related  to  the  windspeed  than  was  the  algebraic  value. 
Overall,  there  was  a  tendency  for  a  better  correlation  between  the  wind- 
speed  and  acid  at  Intermediate  radii  (viz,  3°,  4°,  and  5°).  However,  most 
of  the  correlations  at  separate  radii  were  less  than  0.2.  The  correlation 
between  Tgg  and  windspeed  was  better  at  the  inner  radii  of  the  storm 
(viz,  R  <  5°). 

The  assumption  that  the  orientation  of  the  cloud  bands  would  be  closely 
correlated  with  wind  direction  was  tested  by  examining  the  relation  of  ac^^ 
to  a0bg  and  a  i .  The  highest  correlations  were  generally  found  at 
R '  _>  3  ,  where  the  correlation  coefficients  at  separate  radii  were  about 
0.4  for  a0ks  and  0.2  for  are^.  Errors  in  the  wind  direction  may  have 
contributed  to  the  poor  correlations. 

The  a0bs  and  are^  were  also  correlated  with  windspeed,  the  coefficients 
being  less  than  0.2,  The  windspeed  had  a  slightly  better  correlation  with 
acid  than  with  a0fcs. 

Results  from  the  Multiple  Screening  Regression  Method  showed  that  about 
63%  of  the  variance  of  the  windspeed  could  be  accounted  for  by  the  R~0«5, 
Only  an  additional  2%  to  3%  could  be  accounted  for  by  the  various  forms  of 
TgB  and  ac^^,  which  is  not  considered  significant. 

When  nine  forms  of  the  independent  variables  were  screened  against  the 
variation  of  windspeed  in  the  tangential  direction,  only  about  3%  of  the 
variance  of  AV/A0  was  accounted  for  by  dci^   and  Aa/AR. 

V.   SUGGESTIONS  FOR  FUTURE  RESEARCH 

Further  exploration  of  the  relation  Rx,  such  as  that  studied  by  Arnold 
(1972),  may  be  useful.   In  the  absence  of  wind  observations,  a  simple  form 
of  the  radial  relation  could  be  used  as  a  first  guess  of  the  windspeed 
distribution. 

Eq.  (10)  relates  the  tangential  component  of  the  windspeed  to  the  radius 
and  to  the  wind  direction  crossing  angle,  and  should  therefore  be  tested 
using  only  the  tangential  component. 

Changes  of' the  windspeed  in  time  may  be  better  related  to  the  two 
satellite  parameters  than  is  the  windspeed  itself.   This  is  supported  by 
recent  unpublished  work  of  Ushijima  (1975)  in  which  he  argues  that  the 
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divergence  field  is  one  of  the  most  Important  factors  in  the  formation  of 
the  rainbands.   Because  the  windspeed  profiles  in  the  radial  direction  can 
be  generally  expressed  by  vR^*"  =  constant,  the  radar-observed  rainband 
crossing  angles  should  become  smaller  when  acted  upon  by  such  a  profile. 
Based  on  these  arguments,  the  1-  to  24-hour  changes  of  the  windspeed 
should  be  examined. 

There  are  numerous  ways  in  which  satellite  data  can  be  studied  with 
regard  to  the  windspeed  distribution  problem.   Short-term  changes  of  the 
windspeed  derived  from  cloud  motion  beyond  the  central  dense  overcast  area 
of  tropical  cyclones  can  be  obtained  from  SMS  data  and  should  be  used  along 
with  conventional  data.   It  has  been  suggested  by  Gruber  (1975)  that  the 
perturbation  of  the  infrared  temperature  at  individual  radii  be  tested  for 
possible  relation  to  the  asymmetries  of  the  windspeed.   Other  forms  of  the 
infrared  temperature  field  could  also  be  studied. 

Finally,  studies  similar  to  those  of  Gentry  (1964)  and  Senn  and  Hiser 
(1959) ,  but  with  accompanying  satellite  data,  are  needed  to  shed  more  light 
on  relations  between  the  low-level  convective  cloud  bands  or  radar-observed 
rainbands  and  the  windspeed  in  tropical  cyclones. 
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VIII,   APPENDIX 

DESCRIPTION  OF  SCREENING  REGRESSION  METHODS 

Screening  Regression  methods  are  useful  techniques  in  Isolating  from  a 
set  the  important  or  significant  independent  variables  that  can  be  used 
as  predictors.   In  this  study,  first  a  Stagewise  Screening  Regression 
Method^  was  utilized  (Courtesy:  F,  W.  Nagle,  NESS);  subsequently,  a  Stepwise 
Multiple  Regression  Method2  was  applied  (Courtesy:   C.  M.  Hayden,  NESS). 

The  first  step  of  both  methods  is  to  calculate  the  individual  simple  linear 
correlation  coefficients  (r^n)  between  the  dependent  variable  X^  and  each  of 
the  independent  variables  X2,X3, . . . . ,Xn.   The  variable  X^  is  first  selected 
such v  that   2  o 

r_,  >  r-jjj  from  this  point,  the  two  procedures  differ. 

In  the  Stagewise  method,  the  error  (E^^)  about  the  regression  line 
Xi  =  bik^k  is  then  treated  as  the  dependent  variable  with  the  foregoing 
process  repeated  until  there  is  no  significant  reduction  in  the  explained 
variance  of  the  dependent  variable  (E^#j[). 

2 

In  the  Stepwise  Multiple  Regression  method,  once  the  largest  r^  has  been 

selected,  the  partial  correlation  coefficients,  r^.k  for  i,k  -  2 n 

and  I  /  k  are  calculated.  The  variable  Xg  is  then  selected  such  that 

2       2 
r,  .  >  r^  k  where  g,k  =  2 n  and  g  |*  i,k.  The  process  of  selecting 

the  variable  with  the  highest  partial  correlation  coefficient  is  continued 
until  an  "F  test"  shows  no  significance  or  all  of  the  independent  variables 
have  been  selected. 

Although  the  Stagewise  method  allows  for  a  repetitive  selection  of  the 
same  independent  variable  on  alternating  passes,  while  the  Stepwise 
Multiple  Regression  process  allows  for  the  selection  of  each  independent 
variable  only  once,  the  results  obtained  by  applying  both  methods  were  very 
similar.  Certainly,  the  Stagewise  method  is  a  logical  procedure,  but 
because  it  operates  directly  on  the  residual  error  (E^^)  after  the  first 
pass,  it  is  not  a  statistically  rigorous  method,  as  Is  the  Stepwise  Multiple 
method,  which  operates  directly  on  the  variance  explained  by  each  of  the 
independent  variables. 
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